Enumeration of States in a Periodic Glass 
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We present an analytic enumeration of the metastable states, N s , in a periodic long-range Joseph- 
son array frustrated by a transverse field. We find that the conngurationai entropy, S con f = In N s , is 
extensive and scales with frustration, confirming that the non-random system is glassy. We also find 
that S con f is different from that of its disordered analogue, despite that fact that the two models 
share the same dynamical equations. 



The problem of vitrification has been studied for many 
years and remains a topic of much current interest |Q. 
Recently there has been significant progress in our con- 
ceptual understanding of glass formation in the absence 
of quenched disorder due to studies of exactly soluble 
models . Many features of these regular microscopic 
systems are the same as those of known disordered Hamil- 
tonians || . More specifically, the dynamical equations of 
these periodic glasses are identical to those of (intrinsi- 
cally random) spherical spin glass models; many signa- 
tory properties such as history-dependence and ageing 
follow from these equations. These results suggest that 
one can study periodic glasses via a mapping to disor- 
dered ones for which analytical tools are well-developed 
P,^0[. Here, however, we show that some physical prop- 
erties of periodic and disordered models which share the 
same dynamical equations are different. In particular 
we have calculated the number of metastable states, N s , 
with different ordering of the superconducting phases in 
a periodic long-range Josephson array. We find that the 
configurational entropy, S con f = In N s , is extensive prov- 
ing that this system is indeed in a glassy phase. However, 
Sconf in this model is distinct from the configurational 
entropy of its disordered counterpart, despite the fact 
that they share the same dynamical equations. 

We shall consider the following periodic model, which 
has the advantage that it can be studied both analyt- 
ically HJll[l and experimentally The proposed 
array is a stack of two mutually perpendicular sets of 
N parallel wires with Josephson junctions at each node 
that is placed in an external tranverse field, H. The 
classical thermodynamic variables of this system are the 
superconducting phases associated with each wire. Here 
we shall assume that the Josephson couplings are suffi- 
ciently small so that the induced fields are negligible in 
comparison with H. We can therefore describe the array 
by the Hamiltonian 
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with J,- 



^-exp(2majk/N) and 1 < (J, k) < N 



where J mn is the coupling matrix 



where j(k) is the index of the horizontal (vertical) wires; 
Zm = e ^ m where 1 < m < 2N and the <f> m are the phases 
associated with the superconducting order parameters of 
the 2N wires. Here we have introduced the flux per unit 
strip, a = NHl 2 /$o, where I is the inter-node spacing 
and $0 is the flux quantum; the normalization has been 
chosen so that Tq does not scale with N. 

Because every horizontal (vertical) wire is linked to 
every vertical (horizontal) wire, the number of nearest 
neighbors in this model is N; we can therefore study it 
with a mean-field approach. Such an analysis of the ther- 
modynamic properties indicates that at sufficiently low 
temperatures the paramagnetic phase becomes unstable 
with aN modes becoming simultaneously degenerate [0 . 
We speculated previously that any linear combination of 
these modes would result in a metastable state at lower 
temperatures leading to S con f ~ aN. The results of the 
explicit calculation described below confirm this conjec- 
ture. 

In the limit of small transverse field < a < 1), the 
dynamical equations of the periodic array are identical to 
those of the p — 4 (disordered) spherical spin glass model 
p| . These equations indicate that all metastable states 
are formed at the transition, with no further subdivision 
occurring at lower temperatures [ pT[ . This conclusion 
agrees with that found from a direct study of the TAP 
solutions of the p — 4 spherical spin glass model . It 
is thus sufficient to calculate the number of states, N s , in 
the periodic model at zero temperature. Because there is 
no average over disorder in this regular array, there is no 
distinction between ln{iV s ) and (lniV s ) as occurs for in- 
trinsically random systems, and hence no corresponding 
problems of relevance or interpretation of ]xi{N s ) that one 
usually calculates instead of the more physical (lniV s ). 

At zero temperature the defining characteristic of a 
metastable state is that each spin should be parallel to 
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its associated molecular field M|. This physical con- 
dition results in a highly nonlinear equation which de- 
termines the number of states. Previous studies of dis- 
ordered xy spin systems indicate that the crucial non- 
linearities appear in the expression for the amplitude, 
due to the quasi-linear nature of the phase component. 
More specifically, the configurational entropy for the xy 
Sherrington-Kirkpatrick model jl?]] is different from that 
of its Ising counterpart by a numerical factor of O(l). 
This observation allows one to reduce the array problem 
to that of Ising spins, s m — ±1, thereby simplifying the 
technical presentation. In this case, the number of states 
is given by 
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where we have used an integral representation of the 9- 
function, V(j> = J] m 2 ti£ m and = (^m + 4>n)- 

We use the cumulant expansion to determine the inte- 
gral (^) . Since only terms with an even power of each spin 
variable contributes to this sum, it can be represented as 
a closed loop product of ( JA)ij matrices. Moreover the 
structure of the periodic array is such that only loops 
containing an even product of these matrices are possi- 
ble ||,|l9|]. For the explicit calculation we shall need the 
moments of the couplings (J 2p )- We know from previous 
work on the array ]qjl9|] that 
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(J 2p ) 
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a p-l N 2p-l ( 5 ) 

where Jo = 1; this result can be understood physi- 
cally since only horizontal (vertical) wires separated by 
d < l/a contribute coherently to this sum independent 
of disorder. We emphasize that so far our calculation 
can be applied to both periodic and random arrays; the 
only difference is that in the latter case the couplings are 
averaged over disorder. 



The summation over spin variables in can be graph- 
ically represented by the sum of diagrams with closed 
loops, shown in Fig. 1. Diagrams with intersecting 
loops (see Fig. 1) are small in a due to the structure 
of the moments (|^); thus we can neglect such diagrams 
for a < 1 which simplifies the problem. Furthermore 
each closed loop diagram, (see Fig. la), makes a con- 
tribution K 2 p = j T>4> Tr {JA) 2p to (|) where p refers 
to the number of nodes in the graph. The integration 
over the phase variables makes all sites equivalent, so 
that K 2 p = (J 2 p) J V<t> Tr A 2p . Summing all the loop 
diagrams and rescaling cf> — > ^/acj), we have 



N S =2 N J V(j)e aNS ^ 



(6) 



More explicitly, because Aij = (<j)i + <j>j ) 
Tr A™ = £ (ftx+fcjfafa (7) 



In each term of the resulting polynomial each index is 
represented zero, one or at most two times; the resulting 
contributions to -^-TrA" are Yli 1 — 1> ~k Si 0» an d 
~k Si 0? respectively. Therefore ^j-TrA 2p is a function 
of only t = -k ^ <f> and u = ^2 ( t >2 ■ I n order to calculate 
TrA™ explicitly, we separate the two terms in which each 
index appears once and only once from those in which 
some indices are repeated and others do not appear at 
all. Furthermore, it is convenient to compute the second 
term via its derivative: 



Tr A" = 2f 



u dTrA r ' 
du 
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Explicitly the derivative in the second term is given by 
<9Tr A 



du 



= nNC n -2 
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FIG. 1. Diagrams for the spin trace: (a) Closed 
non-intersecting loops that give the main contribution 
to N s in the limit a < 1 and (b) Intersecting loops that 
are next order in a. 



Here C n can be determined recursively, once again sep- 
arating into terms which involve only distinct and re- 
peated indices. Explicit computation of the sum in equa- 
tion (H), involving some algebra, yields 
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m{(l+i 2 ) 2 + 2u(l - t 2 ) + u 2 } . (11) 
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Because S depends on <f> solely via t and u denned above, 
the TV-dimensional integral over cf> in (JsJ) can be con- 
veniently calculated by including two additional inte- 
grals over t and u, and the factors S (t — < p) an d 




FIG. 2. S(t, A) for a typical a (a 
the saddle-point. 



; 0.1) which displays 



We may use the integral representation of the 6- 
functions and integrate over each (pi independently; the 
result can be concisely presented using the function 
Q(x) = §{l+Erf(§)} so that the integral in (§) be- 
comes 



IOC pOO f / \ *\ 

\ / ^\ / M A.. dr e N{aS+irt+Au+lne(^)} 



N s = 2 JV / dX dt du 
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where S is given by (|ll|) . The factor of TV in the expo- 
nential of the integral in ( |l2|) allows us to evaluate this 
integral by a saddle-point approximation. We shift our 
contour in t by it to look for a solution, since on physical 
grounds we expect the saddle-point to be real. This is 
indeed the case, as displayed in Figure 2, and the result 
is 



^conf 



IniV, 



Na 1 
~ n 2^' 



(13) 



for a<l. This equation for the configurational entropy 
is the main result of this paper. It indicates that S con f is 
extensive, as expected for a glass. Furthermore S con f is 
proportional to aN, consistent with the fact that there 
is vanishing frustration in the limit of a < Physi- 
cally, we note that the configurational entropy in ( [13] ) is 
proportional to the effective number of spins aN in the 
periodic model; in this array, phases (and hence spins) 
are correlated on a length-scale — where I is the internode 
spacing. Eq. ([l3]) also concurs with a conjecture based 
on a previous stability analysis for the high temperature 
paramagnetic phase 



The calculation above can be generalized to find the 
number of states in which all molecular fields are greater 
than a threshold value, hth- For small hth the asso- 
ciated configurational entropy decreases as 8S con f ~ 
—No^l^hth- We note that both S con f and SS con f are 
proportional to N. Empirically we believe that states 
in which some spins experience a very small molecular 
field are only marginally stable. Thus this perturbative 
change in the configurational entropy induced by finite 
hth ^ 1 implies that in the thermodynamic limit the 
number of marginal states is extensive. 

We can also extend this calculation to determine the 
number of states as a function of their energy. To do 
this we introduce an additional S(NE — J2mn ^mnSmS n ) 
function in the definition (||) which extracts only states 
with energy E. Using an integral representation of this 
5 function with an additional variable fj, we see that T>4> 
in Eq (|J) becomes T>4> = n m 2 -Kiftm+^) an< ^ ^ e ex P onen ~ 
tial acquires an additional contribution iN^tE. Finally in 
the final integral (|l2| ) the exponential acquires the term 
iNEu. and function 0(t/\/A) is replaced by 
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We evaluate the resulting integral by a saddle-point ap- 
proximation. Three out of the five resulting saddle-point 
equations can be solved analytically; the remaining two 
contain an error function and must be solved numerically. 




FIG. 3. Sconf (E) for a = 0.1 

Our numerical solution for the configurational entropy 
as a function of energy is displayed in Figure 3 for a = 
0.1. Qualitatively it is similar to the behavior of S con f{E) 
for spin glasses with intrinsic randomness fl6|-|Ts|]. How- 
ever its analytical structure is distinct, specifically from 
that of the p = 4 (disordered) spherical model p8[ . Fur- 
thermore all states displayed in Fig. 3 are stable, whereas 
in the p-spin model all states above a certain threshhold 
energy are unstable. We note, however, that here we are 
only considering stability with respect to single-spin flips; 
by contrast in the p-spin model stability is defined with 
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respect to continuous deformations. 

The result ( |l3| ) was derived in the limit of small a, and 
we have also checked the qualitative validity of our con- 
clusion for a = 1 with xy spins using a different method 
than that discussed above. We emphasize that the a = 1 
case is very special, since in this instance is a unitary 
matrix. Furthermore we note that for a = 1 the cou- 
pling matrix J^- is identical to a discrete Fourier trans- 
form; then the condition that the Fourier transform is 
flat defines a stable spin configuration. It is convenient 
to write this condition as B r = X^feLo z k z k+r equals zero 
for r 7^ 0. Then we obtain the equation for the number 
of states 



N 

N s = / Y\ dz r dz. 



'- 6(B r ) det— . 15 

Z7T OZ 



Since the B r are pseudo-random variables with (B 2 ) = N 
we may assume that they obey a Gaussian distribution 
with P(B) ~ 2TrN e ~^W^~- ^he determinant in jl5| ) can 
be evaluated exactly by noting that — z* +k imply- 
ing that for z's that satisfy the condition B r = all rows 
are orthogonal. The length of each row is N so that the 
determinant = N N . Therefore the average over B 



in (15) results in 



Kirkpatrick model, even though they share the same dy- 
namics JT^ I ; however S con f for the periodic and the ran- 
dom networks are identical, which may be coincidental. 
These qualitative results remain valid for a = 1, thus sug- 
gesting that they are independent of our approximation 
given the absence of commensurability. We note, how- 
ever, that the dynamical equations change their form at 
a ~ 1, and become similar to those of disordered spher- 
ical spin glass models with interactions that are sums of 
couplings with p > 4; qualitatively this does not seem to 
affect the number of states. 
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S conf = lnN s =N\n- 



(16) 



FIG. 4. Numerical results for the configurational en- 
tropy, s = S con f /N, for a = 1 as a function of N. Note 
that for N < 12 the possible error is less than the size of 
the points. 



This expression was derived in the limit of N ^> 1 ; from 
its form it is clear that this asymptotic result can be at- 
tained only at numerically large N (at least N ^> 2tt) . In 
order to check the behavior of S con f for a = 1 at moder- 
ate N, we have performed a direct numerical minimiza- 
tion of the Hamiltonian in (Q); the results, displayed in 
Fig. 3, clearly indicate that the configurational entropy 
is extensive. We expect that S con f crosses over from 
Sconf ~ Nh\2 (obtained numerically as shown in Fig. 
3) to the analytical result (|l6|) at N w 20; however from 
a practical standpoint this crossover is inaccessible nu- 
merically due to the extensive nature of the ground-state 
manifold. 

To summarize our results, we have calculated the num- 
ber of states in a periodic glass and have found that the 
configurational entropy, Sconf = hiiV s , is proportional 
to aN (Eq. @). It is thus extensive but is different 
from that of the p = 4 (disordered) spherical spin glass 
model JtJ despite the fact that their dynamical equations 
are identical [jn 14|. Furthermore, the complexities (the 
number of states as a function of energy) of these two 
models have different functional forms. It appears that 
rather different microscopic models can have the same 
dynamical behavior, as evidenced by the fact that the 
parameter a (for a < 1) does not affect the rescaled 
response. As an aside, we note that S CO nf for the disor- 
dered array is also distinct from that of the Sherrington- 



The fact that the periodic long-range array and the 
p = 4 spherical (disordered) spin glass have different 
configurational entropies suggests that the detailed struc- 
ture of their respective phase spaces, e.g. barriers, basins 
of attraction, is quite distinct. This is not inconsistent 
with previous work which indicates that these two mod- 
els share the same dynamical equations Jll],[l4| . The lat- 
ter were derived in the thermodynamic limit (JV — > oo), 
and do not describe transitions between metastable states 
that might occur at finite N. Thus these equations only 
probe the system's responses in the vicinity of a typical 
state where it is trapped, and are insensitive to the to- 
tal number of metastable configurations. The distinction 
between the configurational entropies of the periodic and 
disordered models indicates differences in their physical 
properties at finite N and, most likely, in the charac- 
teristics of analogous finite-range systems. We therefore 
believe that a detailed study of the structure of their re- 
spective phase spaces could give insight about the physics 
of finite-range problems and about differences between 
disordered and structural glasses. 

Experimentally states are accessed with some weigh- 
ing factor which is ignored by the quantity N s calculated 
here; for example in thermal equilibrium they would be 
weighted with a Boltzmann factor, whereas in a rapid 
quench they would be weighted according to the size of 
their basins of attraction. Furthermore it seems likely 
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that the distribution of the basins of attraction is 
broad; if so, some states would have weights much larger 
than others in a rapid quench. Nonetheless we believe 
that the main conclusion of this paper remains valid; the 
number of states can be very different for systems which 
share the same dynamical equations. 
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